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METHOD OF COMBINING VERTICAL AND MAGNETIC DIPOLE 
INDUCTION LOGS FOR REDUCED SHOULDER AND BOREHOLE EFFECTS 

BACKGROUND 

Field of the Invention 

[0001] The present invention generally relates to methods and systems for measuring the 

properties of a formation through which a borehole passes. More particularly, the present 
invention relates to induction logging tools that measure the formation response to vertically- 
and horizontally-oriented dipoles. Still more particularly, the present invention relates to a 
method for combining the measured responses to achieve an induction.log having reduced 

'•3 shoulder and borehole effects. 

m 

m 

■■isJ 

:i| Description of the Related Art 

CI?! 

,[0002] The basic principles and techniques for electromagnetic logging for earth formations are 

a 

l II well known. Induction logging to determine the resistivity (or its inverse, conductivity) of earth 
UH formations adjacent a borehole has long been a standard and important technique in the search 
^ for and recovery of subterranean petroleum deposits. In brief, the measurements are made by 
inducing electrical eddy currents to flow in the formations in response to an AC transmitter 
signal, and measuring the appropriate characteristics of a receiver signal generated by the 
formation eddy currents. The formation properties identified by these signals are then recorded at 
the surface as a function of the depth of the tool in the borehole. This record of the measurements 
is commonly termed "a log", or more specifically, "an induction log" of the formation. 
[0003] It is well known that subterranean formations surrounding an earth borehole may be 

anisotropic with regard to the conduction of electrical currents. The phenomenon of electrical 
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anisotropy is generally a consequence of either microscopic or macroscopic geometry, or a 
combination thereof, as follows. 
[0004] In many sedimentary strata, electrical current flows more easily in a direction parallel to 

the bedding planes than it does in a direction perpendicular to the bedding planes. One reason is 
that a great number of mineral crystals possess a flat or elongated shape (e.g., mica or kaolin). At 
the time they were laid down, they naturally took on an orientation parallel to the plane of 
sedimentation. The interstices in the formations are, therefore, generally parallel to the bedding 
plane, and the current is able to easily travel along these interstices which often contain 
electrically conductive mineralized water. Such electrical anisotropy, sometimes called 
^ microscopic anisotropy, is observed mostly in shales. 



different lithological characteristics and, therefore, different resistivities. In well logging 



\$\ systems, the distances between the electrodes or antennas are great enough that the volume 
j;3 involved in a measurement may include several such thin beds. When individual layers are 
: neither delineated nor resolved by a logging tool, the tool responds to the formation as if it were 
j;;f j a macroscopically anisotropic formation. A thinly laminated sand/shale sequence is a particularly 

important example of a macroscopically anisotropic formation. 
[0006] If a sample is cut from a subterranean formation, the resistivity of the sample measured 

with current flowing parallel to the bedding planes is called the transverse or horizontal 
resistivity Ph. The inverse of p H is the horizontal conductivity a H . The resistivity of the sample 
measured with a current flowing perpendicular to the bedding plane is called the longitudinal or 
vertical resistivity, pv, and its inverse the vertical conductivity oy. The anisotropy coefficient X 

is defined as: X = ^G h /o v . 




[0005] Subterranean formations are often made up of 



a series of relatively thin beds having 
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[0007] In situations where the borehole intersects the formation substantially perpendicular to 

the bedding planes, conventional induction and electromagnetic wave propagation well logging 
tools are sensitive exclusively to the horizontal component of the formation resistivity. This is a 
consequence of the induced currents flowing in horizontal planes in the absence of formation dip 
or well deviation. Indeed, regarding Galvanic devices, the lack of sensitivity to anisotropy is 
even more stringent due to the "paradox of anisotropy", which states that any array of electrodes 
or sensors deployed along the axis of a wellbore in a vertical well is insensitive to the vertical 
component of resistivity, despite the intuitive expectation to the contrary. 

[0008] However, it becomes possible to measure the vertical resistivity by orienting antenna 

coils away from the axis of the induction tool. An example of a commercial instrument that 

;5j measures both horizontal and vertical resistivity is described by B. Kriegshauser, et al., describe 

A 

hij this instrument in "A new multicomponent induction logging tool to resolve anisotropic 
ill formations", 41 st Annual Logging Symposium, Society of Professional Well Log Analysts, paper 

O D, pps. 1-14, 2000. This instrument employs multiple multi-component coils (i.e. transmitter and 

\%\ 

: ;*| receiver coils having axial and transverse orientations). Other tools designed to measure both 
{4i horizontal and vertical resistivity are described in U.S. Patent No. 4,302,723 entitled "Apparatus 
and method for determining dip and /or anisotropy of formations surrounding a borehole" by J. 
Moran and in U.S. Patent Application No. 09/583,184, entitled "Method for Iterative 
Determination of Conductivity in Anisotropic Dipping Formations" and filed May 30, 2000, by 
inventors L. Gao and S.C. Gianzero. 
[0009] Unfortunately, induction logging of both horizontal and vertical resistivities suffers from 

what is termed the "shoulder effect". The true formation resistance is believed to vary as a 
function of depth in a fairly rectilinear fashion, i.e. the resistance changes discontinuously as one 
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crosses boundaries between formation layers. However, because the tool measures the properties 
of a nonzero formation volume, the tool measurement actually varies gradually as the tool moves 
across a boundary. This is because some portion of formations on both sides of the boundary are 
within the measurement volume and consequently both contribute to the measurement until the 
tool is well past the boundary. In induction tools the measurement volume extends infinitely in 
all directions, with the effect of the formation resistance falling off exponentially with distance 
from the tool. The contribution of the formation outside the area of immediate interest produces a 
generally undesirable "softening" of the induction log, i.e. the resistance changes gradually as 
boundaries are encountered. This effect is most evident when thin formation layers are 
J;* considered. Where an ideal resistivity tool would show two discontinuities from the opposite 
boundaries of the thin layer, the actual log shows only a small bump between the boundaries. In 

■So 

f n ) effect, the "shoulders" of the ideal resistivity tool have been erased. 

=P010] Also, an induction log of vertical resistivities tends to suffer from an inordinate borehole 

Q effect. That is, the borehole fluid and borehole geometry affect the measurement of the formation 
resistivity in an undesirable way. A method of induction logging that reduces both the shoulder 

i : rj 

effect and the borehole effect would be very desirable. 



SUMMARY OF THE INVENTION 
[0011] Accordingly, there is disclosed herein a method of enhancing the vertical resolution of an 

induction tool, in a manner that may advantageously also reduce undesirable borehole and 
"negative resistivity" effects. In one embodiment, the method comprises: a) obtaining a vertical 
magnetic dipole (VMD) response signal from a transmitter-receiver array of antenna elements 
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having magnetic dipoles oriented parallel to a tool axis; b) obtaining a horizontal magnetic dipole 
(HMD) response signal from a transmitter-receiver array of elements having magnetic dipoles 
oriented perpendicular to the tool axis; and c) combining the VMD and HMD response signals to 
obtain a combination response signal. When the relative weights of the VMD and HMD response 
signals are set as described herein, the combination response signal (and any log calculated 
therefrom) has a narrow, substantially rectilinear, vertical measurement profile. Further, the 
combination response signal is relatively insensitive to the negative borehole contribution 
intrinsic to HMD. 

BRIEF DESCRIPTION OF THE DRAWINGS 
012] A better understanding of the present invention can be obtained when the following 

detailed description of the preferred embodiment is considered in conjunction with the following 
drawings, in which: 

Fig. 1 shows a typical logging-while-drilling (LWD) environment; 
Fig. 2 shows a coil configuration of a conceptual triaxial induction tool; 
Fig. 3 is a graph of a vertical measurement profile for a one transmitter, two receiver, 
HMD array; 

Fig. 4 is a graph of a vertical measurement profile for a one transmitter, two receiver, 
VMD array; 

Fig. 5 is a graph comparing a vertical measurement profile of a combined HMD/VMD 
tool to that of a VMD-only tool and HMD-only tool; 

Fig. 6 is a graph of an integrated radial measurement profile for a one transmitter, two 
receiver HMD array; 



64536.03/1391.30200 



-5- 



Fig. 7 is a . graph of an integrated radial measurement profile for a one transmitter, two 
receiver VMD array; 

Fig. 8 is a graph of an integrated radial measurement profile for a combined HMD/VMD 
tool; and 

Fig. 9 is a graph showing the sensitivity of the combined HMD/VMD tool to formation 
anisotropy. 

♦ 

[0013] While the invention is susceptible to various modifications and alternative forms, specific 

embodiments thereof are shown by way of example in the drawings and will herein be described 
in detail. It should be understood, however, that the drawings and detailed description thereto are 
not intended to limit the invention to the particular form disclosed, but on the contrary, the 

m 

\/\ intention is to cover all modifications, equivalents and alternatives falling within the spirit and 

:** 

!)H scope of the present invention as defined by the appended claims. 

a 

ill 

O DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

m 

Tool Configuration 

[0014] Turning now to the figures, Fig. 1 shows a well during drilling operations. A drilling 

platform 2 is equipped with a derrick 4 that supports a hoist 6. Drilling of oil and gas wells is 
carried out with a string of drill pipes connected together by "tool" joints 7 so as to form a drill 
string 8. The hoist 6 suspends a kelly 10 that lowers the drill string 8 through rotary table 12. 
Connected to the lower end of the drill string 8 is a drill bit 14. The bit 14 is rotated (and drilling 
accomplished) by rotating the drill string 8, by use of a downhole motor near the drill bit, or by 
both methods. 
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[0015] Drilling fluid, termed mud, is pumped by mud recirculation equipment 16 through supply 

pipe 18, through drilling kelly 10, and down through the drill string 8 at high pressures and 
volumes to emerge through nozzles or jets in the drill bit 14. The mud then travels back up the 
hole via the annulus formed between the exterior of the drill string 8 and the borehole wall 20, 
through a blowout preventer (not specifically shown), and into a mud pit 24 on the surface. On 
the surface, the drilling mud is cleaned and then recirculated by recirculation equipment 16. The 
drilling mud is used to cool the drill bit 14, to carry cuttings from the base of the bore to the 
surface, and to balance the hydrostatic pressure in the rock formations. 

[0016] For LWD, downhole sensors 26 are located in the drill string 8 near the drill bit 14. The 

sensors 26 preferably include an induction tool having multi-axial transmitters and receivers. In a 

M 

r'ffj preferred embodiment, downhole sensors 26 are coupled to a telemetry transmitter 28 that 

Si 

gij transmits telemetry signals by modulating the mud flow in drill string 8. A telemetry receiver 30 

?X\ is coupled to the kelly 10 to receive transmitted telemetry signals. Other telemetry transmission 

i|3 techniques are known and may be used. The receiver 30 communicates the telemetry to a surface 

? y 

: : 2 installation (not specifically shown) that processes and stores the measurements. The surface 
;|| installation typically includes a computer system of some kind, e.g. a desktop computer. 
[0017] The drill bit 14 is shown penetrating a formation having a series of layered beds 32 

dipping at an angle. A first (x,y,z) coordinate system associated with the sensors 26 is shown, 
and a second coordinate system (x",y",z") associated with the beds 32 is shown. The bed 
coordinate system has the z" axis perpendicular to the bedding plane, has the y" axis in a 
horizontal plane, and has the x" axis pointing "downhill" along the bedding plane. The two 
coordinate systems are related by two rotations, known as the strike angle and the dip angle. Any 
vector in one of the coordinate systems can be expressed in the other coordinate system using 
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basic rotational transform matrices. Consequently, given measurements in the coordinate system 
of the induction tool, the corresponding measurements in the coordinate system of the beds can 
be determined if the dip and strike angles are known. Methods for doing this are discussed in 
U.S. Patent Application No. 09/925,997, "Virtual Steering of Induction Tool for Determination of 
Dip Angle" and filed August 9, 2001, by L. Gao and S.C. Gianzero. 
[0018] In the preferred embodiment, the induction tool employs multi-axial transmitter and 

receiver "triads". Transmitter-receiver coupling measurements may be made along each axis and 
between axes as well. Because the principle of linear superposition applies to electromagnetic 
fields, rotational transforms can be used to manipulate the coupling measurements. The 

\* measurements of "virtual" transmitters and receivers having arbitrary orientations can be 

Q 

2 synthesized in this manner. So, despite the actual orientation of the tool transmitters and 

i.y receivers with respect to the formation bedding plane, it is possible to determine the 

iljfl measurements that would be made by transmitters and receivers that line up with the coordinate 

||3 system of the formation layers. To simplify the following discussion, it is assumed that indeed, 

* iJ 

such measurements are being made by the tool. 
i|j)019j Fig. 2 shows a conceptual sketch of a coil arrangement for a downhole induction tool. A 

triad of transmitter coils T x , T y and T z , each oriented along a respective axis, is provided. Two 
triads of similarly oriented receiver coils (R Xi ,R yi , R Z] ) and (R X2 > R y 2 > R z 2 ) are also 
provided, separated from the transmitter triad by LI and L2, respectively. Each of the coils in the 
triads is parallel to the corresponding coils of the other triads, and the triads are spaced apart in 
the z-axis direction. 

[0020] One way to study this arrangement of coils is to consider the horizontal and vertical coils 

separately. Transmitter coil T y has an axis that is perpendicular to the axis of the tool. When 
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driven with an AC signal, this coil is equivalent to a magnetic dipole oriented along the axis of 
the coil, so this coil is termed a "horizontal magnetic dipole" transmitter. Receiver coils R y . 

respond to magnetic fields in the same manner as magnetic dipoles oriented along the axes of the 
coils, and are accordingly termed "horizontal magnetic dipole" receivers. Transmitter coil T z , 
and receiver coils R z . , are termed vertical magnetic dipole transmitter/receivers because their 

axes are parallel to the axis of the tool. 
Geometric Factors and Vertical Measurement Profile 
[0021] When an AC signal is provided to a transmitter coil T y , the receiver coils R y . provide an 

y t output AC signal in response. For any given signal frequency, the output signal is proportional to 
i;3 the transmitter signal (once a time delay is taken into account). The proportionality constant 

m 

i,f r \ depends on the magnetic coupling between the transmitter coil and receiver coils, which in turn 
■;P is affected by the properties of the formation enclosed by the measurement volume. As 
mentioned previously, the effect that any given portion of the formation has on the measurement 

'•■zj 

'A is dependent on its position relative to the coils. The proportionality constant can be calculated in 

m 

H the following manner: 

ry 

P = a rSt^Sd<pdrdz 

where a is a constant, a is the conductivity of the formation (and may vary within the formation), 
and g is a geometric weighting factor that accounts for the relative contribution of each portion 
of the measurement volume to the overall measurement. The preferred coordinate system here is 
cylindrical, with <p being the azimuthal angle about the tool axis, r being the radial distance from 
the tool axis, and z being distance along the tool axis. The geometric factor for a horizontal 
magnetic dipole transmitter-receiver pair is: 
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Jz T -z R \ r[r 2 +2(z-z r )(z-z R )] 

Shmd i 3 / 3 

2 Vr 2 +(z-z r ) 2 Vr 2 +(z-z r ) 2 
where the transmitter is located on the tool axis at z T and the receiver is located on the tool axis 
atz R . 

[0022] Typically, the sensitivity of the tool to the formation signal is enhanced by combining the 

measurement of two receivers at different positions along the tool axis. This may be done by 
subtracting the response of R y2 from the response of R yi in a manner that cancels any mutual 

inductance signal from the transmitter. Many induction tools also include an additional 
transmitter positioned symmetrically with respect to the receivers for temperature compensation. 
?3 The geometric factor of one or more receivers to one or more transmitters can be determined by 

'•'rjti 

v f} simply summing the geometric factors (with the appropriate sign) for each pairwise combination 
: ?J of transmitters and receivers: 

G hmd =ZZ Shmd 

a T R 

m 

q For the following discussion, let Ghmd be the total geometric factor for the difference signal of 

m 

i!3 two horizontal magnetic dipole receivers in response to a horizontal magnetic dipole transmitter. 

m 

[0023] The vertical resolution of a tool having a given geometric factor can be determined from 

the vertical measurement profile of the tool. The vertical measurement profile can be calculated 
by integrating the geometric factor over the azimuthal angle <p and radius r. 

Vhmd(z) = $Z n G hmd d<pdr 
Fig. 3 shows this vertical measurement profile. 
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[0024] A similar calculation can be done to determine the vertical measurement profile of 

receiver coils R z . to transmitter T z . . The geometric factor of a single vertical magnetic dipole 
transmitter-receiver pair is: 



Svmd 



\z T -z R \ r 3 



-3 ' 



2 7r 2 +(z-z 7 ) 2 V' 2 + (z-2 fl ) 2 
and when the receiver responses are combined, the total geometric factor is: 

Gymd ~ XX Svmd 
T R 

For the moment, assume that G vm d represents the total geometric factor for the difference signal 
of two vertical magnetic dipole receivers in response to a vertical magnetic dipole transmitter. 
Q The vertical measurement profile of the array is: 

m 

l $ This vertical measurement profile is graphed in Fig. 4. 



Enhancing Vertical Resolution 



I j|)025] The vertical measurement profiles shown in Figs. 3 and 4 both suffer from long "tails" 

Q 

ry that gradually fall off in both directions along the z-axis. These tails allow portions of the 
formation outside the immediate area of interest to interfere with the measurement of properties 
inside the immediate area of interest. Further, the negative portions of the horizontal 
measurement profile make interpretation of those logs a matter of some difficulty. Importantly, 
however, these logs can be combined to produce a vertical measurement profile without tails. 
When the logs are combined in the following manner, the effective vertical measurement profile 
of the tool is that shown in Fig. 5. 

Final Response = ax VMD response -bx HMD response 
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where VMD response is the voltage signal of vertical magnetic dipole receiver array divided by 
the voltage signal provided to the vertical magnetic dipole transmitter array. Similarly, HMD 
response is the voltage signal from the horizontal magnetic dipole receiver array divided by the 
voltage signal provided to the horizontal magnetic dipole transmitter array. The parameters a and 
b are optimization parameters whose values are chosen so as to best enhance the vertical 
characteristics of the final response. For example, in Fig. 5, a = 3/2, b = 1/2. Other values of a 
and b may be used to optimize different aspects of the combination process. 

[0026] The dotted line in Fig. 5 is the vertical measurement profile of the horizontal magnetic 

dipole array, the dashed line is the vertical measurement profile of the vertical magnetic dipole 

"Z\ array, and the solid line is the vertical measurement profile of the horizontal/vertical 

a 

jfj combination. Note that the vertical resolution of the combined log is not subject to the 



jj undesirable "tails". While the discussion here has focused on a single-transmitter, two-receiver 
f .?l tool (with both horizontal and vertical components for each transducer), the above combining 



equation works to eliminate measurement profile tails for any combination of transmitters and 
j'jjj receivers having both horizontal and vertical component measurements. 

; |p027] Of course, a designer may wish to make trade-offs to improve the overall utility of a tool. 

In such a scenario, it may be desirable to combine the horizontal and vertical component logs in 
different proportions to reduce, but not necessarily eliminate, the tails of the vertical 
measurement profile. Such a configuration might offer reduced sensitivity to other undesirable 
effects (e.g. borehole effect) or enhanced sensitivity to desirable effects (e.g. anisotropy). 
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Radial Measurement Profile and Investigation Depth 
[0028] The radial measurement profile of an induction tool can be calculated from the total 

geometric factor. For horizontal and vertical magnetic dipole arrays, the radial measurement 
profile is: 

Rhmdir) = Co Z n G hmd d <P dz 

R vmd( r ) = Co t K G vmd d(pdz 
Although the radial measurement profile itself is informative, it is generally considered more 
informative to examine the integrated radial measurement profile. The integrated radial 
M> measurement profile shows the cumulative contribution of the enclosed measurement volume to 

9 

O the overall measurement, and it converges to one (100%) at large radii: 

m 

U! iRhmd(r)= )o R hmd(p) d P 

m 

IR vmd<< r )=\ ) R vmd(p) d P 

9 

! ; |j0029] Fig. 6 shows the integrated radial measurement profile for the horizontal magnetic dipole 

hf 

ZJ one-transmitter- two-receiver configuration. Note the negative region between r=0 and r=l. This 

ry 

negative contribution is undesirable and it causes a significant borehole effect. The integrated 
radial measurement profile for the vertical magnetic dipole configuration is shown in Fig. 7, and 
this curve has a more desirable shape. Unlike the vertical measurement profile case, a long tail is 
considered beneficial for the horizontal measurement profile. Fig. 8 shows the integrated radial 
measurement profile for the combined horizontal and vertical magnetic dipole log. Note that the 
negative contribution from the borehole and region immediately beyond are completely absent. 
[0030] The investigation depth of a tool is defined to be the depth at which the integrated radial 

measurement profile equals 0.5. The investigation depth of the combined horizontal/vertical log 
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is less than those of the vertical magnetic dipoles, but this may be a small price for the greatly 
enhanced vertical resolution of this tool. This is particularly true since the investigation depth for 
the horizontal magnetic dipole is deceptively deep due to the negative portion of the radial 
measurement profile. 
Anisotropy 

[0031] It is known that the horizontal magnetic diople transmitter-receiver pairs are sensitive to 

resistivity anisotropy in the formation. (See, e.g., J. H. Moran and S. Gianzero, Geophysics, Vol. 
44, P. 1266, 1979.) Although some sensitivity is lost, the combination log still retains adequate 
sensitivity to anisotropy. Fig. 9 shows apparent conductivity as calculated from the combination 
!•» log measurements (combo c a ) minus apparent conductivity as calculated from vertical magnetic 

| dipole sonde measurements (VMD c a ). This difference is shown as a function of VMD c a for 

k 

\ : %\ seven different anisotropy coefficients. The spread in the curves indicates the sensitivity of the 
W tool to anisotropy. 

;;fb032] Ultimately, the results of processing the combination and VMD sondes to enhance 

iy 

I'jjj resolution and correct for skin and borehole effects should yield measures of the formation 
conductivities that are closer to the true conductivity values. In effect, their measured 
conductivities should be close to the homogeneous value of formation conductivity. It may be 
desirable to perform such processing first, and then performing the anisotropy determination 
using the more correct values of formation conductivity. 
Conclusion 

[0033] The vertical measurement profile of the combination sonde is compared in Fig. 5 against 

that of the conventional VMD sonde. Clearly the combination sonde is better focused compared 
to both the conventional VMD sonde and the HMD sonde. The large tail outside the transmitter- 
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receiver spacing has been totally eliminated. The negative contribution from the HMD sonde has 



also been removed. Further, the rectangular nature of the combination sonde's vertical 



measurement profile lends itself to deconvolution for improved vertical resolution in the 



formation log. Although skin effect has not been considered here, it is expected that the 



elimination of the shoulder effect will substantially reduce the skin effect, which is itself a far 



field phenomenon. Finally, the combination sonde advantageously retains its sensitivity to 



anisotropy so that anisotropy measurements can be made. 



[0034] Numerous variations and modifications will become apparent to those skilled in the art 

z * 

once the above disclosure is fully appreciated. For example, multiple horizontal magnetic dipole 



!' * array response signals and multiple vertical dipole array response signals may be combined with 
:=1 a total proportionality of vertical to horizontal contributions approximately equal to minus three. 



; y The use of multiple signals may provide improved performance. It is intended that the following 
i!jf! claims be interpreted to embrace all such variations and modifications. 



a 
m 
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